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Modified 50(3) Lattice Gauge Theory at T ^ with Parallel Tempering 
Monopole and Vortex Condensation 
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The deconfinement transition is studied close to the continuum limit of SO(3) lattice gauge theory. 
High barriers for tunnelling among different twist sectors causing loss of ergodicity for local update 
algorithms are circumvented by means of parallel tempering. We compute monopole and center 
vortex free energies both within the confining phase and through the deconfinement transition. We 
discuss in detail the general problem of defining order parameters for adjoint actions. 
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I. INTRODUCTION 



Understanding confinement in SU (N) Yang-Mills the- 
ories remains one of the major challenges of contempo- 
rary particle physics. Lattice simulations have offered 
unique insight into the non-perturbative regularization 
of pure gauge actions transformingunder the fundamen- 
tal representation of SU(N) [J 12], equivalent to the 
quenched limit of full QCD: at non-zero temperature 
they have been shown to possess a phase transition linked 
to the spontaneous breaking of center symmetry [3, Q • 
For N = 2 it is of second order, therefore lying in the 
universality class of the 3-d Ising model. However the 
question whether and in what sense this also holds for 
discretizations transforming under the natural contin- 
uum pure Yang-Mills gauge symmetry group SU(N) /Zjy, 
for N — 2 equivalent to 50(3), still needs to be ap- 
propriately answered According to universality @, 
i.e. expecting the different formulations to be equiva- 
lent in the continuum limit, they should lead to the same 
non-perturbative physics. A discretization which does 
not break the SU(N)/Zn invariance has moreover the 
appeal to preserve the topological properties related to 
■kx{SU{N)/Z n ) = Z N discussed e.g. in 0,1,1). 

Since the lattice link variables gauge transform at dif- 
ferent points U^{x) — > g^(x)Ufj,(x)g(x + fi), SU{N)/Zn 
invariance cannot be recovered from the local cancella- 
tion of the Zjv dependence in g(x) as in the continuum 
and must be imposed directly on U^ix). As a conse- 
quence in adjoint theories regularized on the lattice it is 
by construction impossible to define observables trans- 
forming under the fundamental representation, i.e. sen- 
sitive to the Zjv center of the gauge group: their expec- 
tation value will vanish identically irrespective of the dy- 
namics of the theory. Therefore the symmetry breaking 
arguments for the deconfinement transition mentioned 
above cannot apply. It remains an open question whether 
a non-perturbative regularization of Yang-Mills theories 
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allowing both SU{N) /Zm invariance and non-vanishing 
fundamental observables can be defined. 

In spite of all these interesting problems adjoint actions 
have not been intensively studied in the literature. For 
N = 2 difficulties in their analysis have been well known 
for a long time [IE HH [H, 03 : the theory exhibits a bulk 
transition related to the condensation of Z2 monopole 
charges a c = — 1 which hinders the study of its finite 
temperature properties. First concrete efforts to study 
the theory at finite temperature by implementing a Z2 
m onop ole suppressing chemical potential, as suggested 
in [T2 [TU , were made ten years ago 0, [l5[ reviving the 
interest in the subject. However, given the absence of 
a natural order parameter, attempts to locate a transi- 
tion within phase II (the phase characterized by strong 
Z2 monopole suppression, see Fig. [1]) through thermo- 
dynamic observables were only conclusive in the strong 
coupling region [IE QjJ • In these works it was also first 
observed how in some regions of phase II close to the bulk 
transition the theory possesses new states where the ad- 
joint Polyakov loop La — > —1/3, additionally to the ex- 
pected states where La — > 1. In Ref. 0] a dynamical 
observable measuring the twist expectation value z, i.e. 
the topological index linked to iti(SU(N)/Zn), was in- 
troduced noting that the 5(a c — 1) constraint effectively 
implemented by a Z2 monopole suppression should al- 
low the 50(3) partition function to be rewritten as the 
sum of 5(7(2) partition functions with all possible twisted 
boundary conditions Z\, = \, (i = 0, . . . , 3 for SU(2) on a 
3+1 dimensional torus) [3 M, EE E3 . The L A -> -1/3 
state was thus linked directly to a non-trivial twist expec- 
tation value, equivalent to the creation of a vortex. Cre- 
ating such 't Hooft loop amounts to changing the signs of 
some fundamental plaqucttcs, which however leaves the 
adjoint action unmodified. This implies that AU — in 
the free energy change AF = AU — TAS , which will then 
only receive an entropy contribution. Defining thus the 
't Hooft vortex free energy F/T = - log {Z\ z - X / Z\ Z=Q ) 
simply by the ratio of the partition function in the non- 
trivial twist sector to that in the trivial one, their relative 
weight being measured through an ergodic simulation, 
the 50(3) theory was proposed as the ideal test case to 
check whether the 't Hooft vortex confinement criterion 
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0, [1] could compensate for the absence of an explicit or- 
der parameter linked to center symmetry breaking: in 
the thermodynamic limit (V — iV 3 — > oc) F should van- 
ish in the confined phase while diverging with an area law 
F ~ <jN^ above the deconfinement transition. Working 
without the monopole suppression term proved however 
to be a hurdle, since the "freezing" of twist sectors above 
the bulk transition creates high potential barriers hard 
to overcome even with a multi-canonical algorithm Q, 
making ergodic simulations on top of the bulk transition 
unviable already for volumes larger than 8 3 x 4. Further- 
more, since one would eventually need to go well beyond 
the bulk transition deeply into phase II with the sim- 
ulations, the suitability of multihistogram [22I ] or mul- 
ticanonical methods [23j remains dubious. Ergodicity 
problems and non-trivial twist sectors were not consid- 
ered in Ref. pj,[l7|- 

A particular observation has proven crucial in our tam- 
ing of the tunnelling problem: the bulk transition weak- 
ens with increasing Z2 monopole suppression, eventu- 
ally becoming 2 nd order at some intermediate point [l6j |. 
Through the twist susceptibility the 2 nd order branch of 
the bulk transition was shown to be consistent with the 
4-d Ising model universality clas s fl6l.l24l [2f| [26[ , as ex- 
pected by theoretical arguments |l2lll3j. To actually pin 
down the point where the transition changes from weak 
first to second order is a difficult numerical task. This 
however has no practical consequences, since for the fol- 
lowing it is immaterial whether one deals with a second 
or a very weak first order bulk. 

Although tunnelling among topological sectors is still 
suppressed with a local update algorithm, twists were 
shown to be well defined throughout phase II. La on the 
other hand approximately satisfies (B,[i3| a Haar- measure 
distribution for low Pa, departing from it above some 
P A - The critical value (5 A was seen to scale properly 
with the lattice extent in the Euclidean time direction 
N T 0, [2^| . This hints at a transition line (the dashed 
horizontal line in Fig.[l| separating a confining from a de- 
confining phase in each fixed twist sector [26U281 129I [3(3] 
collapsing on the bulk transition for the N T commonly 
used in simulations. It is therefore sound to conjecture 
that the whole physically relevant 50(3) dynamics lies 
in phase II, the finite temperature transition eventually 
decoupling from the bulk transition for high enough N T 
even without monopole suppression term. Unfortunately, 
according to estimates in Ref. [9( this should not hap- 
pen for lattice sizes smaller than ~ 800 3 x 400. A non- 
vanishing 2,2 monopole chemical potential together with 
an ergodic algorithm suitable for simulations throughout 
phase II seems therefore the only feasible way to gain 
access to the properties of the continuum limit of 50(3). 

Given the failure of center symmetry breaking criteria 
to identify the deconfinement transition in the adjoint 
theory, in [26|, [3l[ the use of the Pisa disorder parameter 
for monopole condensation was proposed. Lines of sec- 
ond order transition properly scaling with N T and end- 
ing on the bulk transition line where actually found at 



fixed twist, with critical exponents consistent with the 3-d 
Ising model. Whether this is the case also for the ergodic 
theory, i.e. summed over all twist sectors, is the subject 
of the present paper. We will employ parallel tempering 
and utilize the mentioned weakening of the bulk transi- 
tion to overcome the high potential barriers that prevent 
tunnelling with local update algorithms. Moreover, er- 
godicity being an essential prerequisite for an unbiased 
measure of the vortex free energy, it is an interesting 
question in its own right whether such observable could 
indeed also play the role of an order parameter for the de- 
confinement transition in 50(3) [t|[32j]. We will extend 
here the discussion of the vortex free energy we have re- 
cently published in Ref. [33| . Some preliminary results 
of the present project were also presented in [341 ]. 



II. THE SETUP AND THE OBSERVABLES 

As anticipated, we will concentrate on the adjoint 
SU(2) Wilson action modified by a Z2 monopole sup- 
pression term 

S = f3 A Y,( 1 -l T ^Up)+\J2(l-°c), (1) 



where U p denotes the standard plaquette variable and 
Tt a O = (Tr F 0) 2 - 1 = Tr F (0 2 ) + 1 the adjoint trace. 
The center blind product a c — Yipedc ^gn-C^FUp) taken 
around elementary 3-cubcs c defines the Z2 magnetic 
charge. Its density M = 1 — (j^ c it c ) tends to unity 
in the strong coupling region (phase I) and to zero in 
the weak coupling limit (phase II), N c denoting the total 
number of elementary 3-cubes. The corresponding path- 
integral quantized lattice theory with the action (jTJ) is 
center-blind in the entire Pa — A plane 12 811 . 

The Pisa disorder parameter \i [35l l36l . l37l . [38j has 
been introduced for action (TTJ in Ref. [31| . Its expecta- 
tion value is given by (fi) — (e~ AS ) , where AS — S M — S 
is the difference of the standard plaquette action 5 and 
an action S M modified by the introduction of an adjoint 
bosonic field transforming at the space boundary under 
G ~ SU(2)/U(X) [35|,[3!|. Its evaluation does not require 
any gauge fixing, a point of view we adopt in what follows 
[33 ]. We want to stress here that the introduction of C* 
boundary conditions in the temporal direction, necessary 
to conserve magnetic charge when defining S M at finite 
temperature, poses no conceptual problem in the adjoint 
theory, being equivalent up to a gauge rotation to a par- 
tial twist, i.e. only in the time direction [3l| . Since our 
adjoint action with periodic boundary conditions allows 
all twist matrices to be generated dynamically in any di- 
rection C* boundary conditions will just amount to 
a relabelling of the twist sectors. We will come back to 
this point later on. 

Appropriate twist variables are introduced by @ 
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Since the temporal twists in the various spatial directions 
Ziiii = 1,2,3 are well identified (cither +1 or -1) for 
each configuration in phase II, the partition functions 
restricted to a fixed twist sector are easy to define as 
expectation values of suitable projectors [8|. Explicitly 
we have 
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z/fc being equal to unity if the configuration belongs to 
the fc th sector and vanishing otherwise. 
From Eq. © it follows that 



F = _Tl 0g J?L 



1 



■log 



3<H>) 



(4) 



The factor in the denominator is due to the three equiv- 
alent ways to measure z ij4 = -1 on T 3 x §, rather than 
one as on § 3 x S; in this way F will be normalized to zero 
if 0- and 1-twists are equally probable. This occurs on 
top of the bulk transition and in some sense everywhere 
in phase I, where twist sectors are however badly defined, 
because of the twist variables ([2]) fluctuating around zero. 

We will employ parallel tempering to achieve crgodic- 
ity over different twist sectors when evaluating the ex- 
pectation values of physical observables, e.g. the Pisa 
disorder parameter and the 't Hooft vortex free energy. 
Simulations have been carried out along the paths shown 
in Fig. [T] The motivation for these choices will become 
clear in the following. Our spatial lattice sizes will vary 
between N s — 12 and N s = 24. The time-like extension 
will remain fixed (N T = 4). 



III. PARALLEL TEMPERING 
A. General description 

In tempering methods some parameters of the action 
are made dynamical variables in the simulations, updat- 
ing the system in an enlarged configuration space. This 
allows a detour in parameter space if a high tunnelling 
barrier is present at some parameter value, resulting in 
an improved algorithm. 

In the method of simulated tempering first proposed 
in 40] the inverse temperature is made a dynamical vari- 
able. With such algorithms considerable improvements 
have been obtained when rendering dynamical e.g. the 
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FIG. 1: Paths chosen for main simulations in the A — (3a plane. 



number of degrees of freedom in the Potts-Model [4l| , the 
inverse temperature for spin glass |42| and the monopole 
coupling in U(l) lattice theory [43|, H|- With dynami- 
cal mass of staggered fermions in full QCD [45[ a better 
sampling of the configuration space has been reported. 
However, simulated tempering requires the determina- 
tion of a weight function in the generalized action, and 
an efficient method of estimating it [42|, |43|, is crucial 
for successfully accelerating the simulation. 

A major progress was the proposal of the parallel tem- 
pering method (PT) [46l l47j , in which no weight function 
needs to be determined. This method has allowed great 
improvements for spin glasses [46j . In QCD with dy- 
namical quark mass better sampling has been reported 
for staggered fermions [48[. In simulations of QCD with 
0(a)-improved Wilson fermions 49] no computational 
advantage has been found when making only two (rel- 
atively small) ho pping p arameter values dynamical. In 
subsequent works [50L l51j with more ensembles and stan- 
dard Wilson fermions a considerable increase of the tran- 
sitions between topological sectors has been observed. In 
Ref. [13] these investigations have been extended to a 
detailed comparison with conventional simulations. Un- 
fortunately no gain could be confirmed in that case due 
to the region of parameter space used in which the mech- 
anism of an easier detour was not available. 

In the present application the fact that above the bulk 
phase transition the barriers between the twist sectors 
cannot be overcome at all by conventional algorithm 
makes PT in any case superior. With a chain of param- 
eter points crossing the transition line along the softer 
branch of the bulk transition the idea of an easier detour 
by tempering is ideally realized. This is also reflected by 
the remarkably good efficiency of PT observed. 
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B. Parallel tempering algorithm 

In standard Monte Carlo simulations one deals with 
one parameter set p and generates a sequence of field 
configurations r{s), where s denotes the Monte Carlo 
time. In our case p will include the coupling /3a and the 
chemical potential A. In parallel tempering (PT) [46ll47| 
one updates K field configurations ? n with n — 1 . . . . , K 
in the same run. The characteristic feature is that the 
assignment of the parameter sets pj with j = 1 , . . . , K 
to the field configurations j= n changes in the course of a 
tempered simulation. The global configuration at time 
s will be denoted by -B(s), ^i(s), F2{ s )t--i fk{ s ) where 
the permutation 

B{a)= fM*) Ms) ... njs) ... n K {s)^ (5) 

describes the assignment of the field configurations 
fn,j(s){ s ) to the parameter sets pj. For short this ap- 
proach is called PT with K ensembles. 

The update of the ? n is implemented through a stan- 
dard Metropolis procedure using the parameter sets pj 
as assigned at a given time. The update of B is achieved 
by swapping pairs according to a further Metropolis ac- 
ceptance condition with probability 

P swap (i,j) =min(l,e- AS ) , (6) 

where the variation 

AS = + Sfa, r ni ) + S{Pj,r nj ) 

- S fa , r n . ) - S (pj , r n . ) (7) 

refers to the action S for the parameter set pj and the 
field configurations p n . . The total update of the Monte 
Carlo algorithm, after which its time s increases by one, 
then consists of the updates of all jr m followed by the full 
update of B with a sequence of attempts to swap pairs. 

Detailed balance for the swapping follows from Eq. ([7]) . 
Ergodicity is obtained by updating all p n and by swap- 
ping pairs in such a way that all permutations of Eq. (JSJ) 
can be reached. There remains still the freedom of choos- 
ing the succession of the individual steps. Our choice 
is such that the updates of all r n and that of B alter- 
nate. Our criterion for choosing the succession of swap- 
ping pairs in the update of B has been to minimize the 
average time it takes for the assignment of a field config- 
uration to the parameters to travel from the first to the 
last pair of parameter values. This has led us to swap 
neighboring pairs and to proceed with this along the re- 
spective path in Fig. [TJ 

Observables of interest, associated to a specific set pj, 
will be denoted as 

OiW^O^wW), j = l,...,K. (8) 

As anticipated above, for the success of the method the 
softening of the bulk transition is crucial, since we need 



to "transport" the tunnelling that occurs in phase I and 
on top of the bulk transition into phase II, where twist 
sectors are well defined but frozen. To work at low A, 
i.e. on top of a strong 1 st order bulk, would select too 
high barriers and kill any hope of ergodicity at large vol- 
ume, as experienced in Ref. [j| [32J. Moreover, in that 
parameter range for lattice sizes reasonably available to 
the simulations (i.e. <C 800 3 x 400!) finite volume effects 
would still cover the physical transition 0] . 

Some care is of course necessary also with our method. 
In particular to maintain a sufficient swapping accep- 
tance rate u>, i.e. to avoid the freezing of twist sectors, 
the distance between neighboring couplings must dimin- 
ish with the volume. On the other hand to keep cross- 
correlations under control one does not wish the accep- 
tance rate to be too high. We have chosen to tune the 
parameters for each path and volume at hand so to keep 
the acceptance rate roughly fixed at around to = 12%, a 
value for which we empirically find a good balance be- 
tween auto- and cross-correlations. We also find that 
performing some standard Metropolis overupdate hits on 
the r n before the actual PT update ([7]) is proposed helps 
in diminishing correlations. 

The relatively low value of lo has an intuitive expla- 
nation: for each parameter set pj one wishes to "ex- 
plore" the various twist sectors for a sufficiently long MC 
time before tunnelling. It causes however also a techni- 
cal problem: if the starting configurations are all in the 
same twist sector the ensemble needs a very long time 
before the "disorder" below the bulk spreads to the con- 
figurations further above it. An efficient way out is to 
randomly choose the twist sectors of the elements in the 
start ensemble. 





= 12 


N a 


= 16 


N s = 


= 20(a) 


N s = 


= 20(b) 


2 x 40000 


2 x 40000 


2 x 40000 


2 x 40000 


A 


Pa 


A 


Pa 


A 


Pa 


A 


Pa 


0.78 


0.960 


0.78 


0.960 


0.77 


0.960 


0.77 


0.960 


0.79 


0.960 


0.79 


0.960 


0.78 


0.960 


0.78 


0.960 


0.80 


0.960 


0.80 


0.960 


0.79 


0.960 


0.79 


0.960 


0.80 


0.975 


0.80 


0.970 


0.80 


0.960 


0.80 


0.960 


0.80 


0.990 


0.80 


0.981 


0.80 


0.970 


0.80 


0.970 


0.80 


1.005 


0.80 


0.993 


0.80 


0.980 


0.80 


0.980 


0.80 


1.020 


0.80 


1.006 


0.80 


0.990 


0.80 


0.990 


0.80 


1.035 


0.80 


1.019 


0.80 


1.000 


0.80 


1.005 


0.80 


1.050 


0.80 


1.032 


0.80 


1.010 


0.80 


1.015 


0.80 


1.065 


0.80 


1.045 


0.80 


1.025 


0.80 


1.035 


0.80 


1.080 


0.80 


1.058 


0.80 


1.040 


0.80 


1.045 


0.80 


1.090 


0.80 


1.070 


0.80 


1.050 


0.80 


1.055 



TABLE I: Lattice sizes, statistics and couplings for PT simu- 
lations in Fig. [7] and Fig. [5] right branch. 

Details for the investigated lattice sizes, the chosen pa- 
rameter sets (A, Pa) j, j = 1, 2, . . . ,K and the statistics 
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N s 


= 12 


N s 


= 16 


N s 


= 20 


2 x 40000 


2 x 40000 


2 x 40000 


A 


0A 


A 


0A 


A 


0A 


0.80 


0.860 


0.80 


0.860 


0.80 


0.860 


0.80 


0.870 


0.80 


0.870 


0.80 


0.870 


0.80 


0.875 


0.80 


0.880 


0.80 


0.875 


0.80 


0.880 


0.80 


0.890 


0.80 


0.880 


0.80 


0.885 


0.80 


0.895 


0.80 


0.885 


0.80 


0.890 


0.80 


0.900 


0.80 


0.890 


0.80 


0.900 


0.80 


0.908 


0.80 


0.895 


0.80 


0.905 


0.80 


0.920 


0.80 


0.908 


0.80 


0.908 


0.80 


0.925 


0.80 


0.925 


0.80 


0.910 










0.80 


0.920 










0.80 


0.925 











TABLE II: Parameter sets for PT runs in Fig. [H left branch. 



N a = 


= 12 


N s = 


= 16 


N s = 


= 20 


N s = 


= 24 


100000 


100000 


100000 


100000 


A 


0A 


A 


0A 


A 


0A 


A 


0A 


0.78 


0.95 


0.78 


0.95 


0.785 


0.95 


0.785 


0.95 


0.79 


0.95 


0.79 


0.95 


0.795 


0.95 


0.795 


0.95 


0.795 


0.95 


0.795 


0.95 


0.7975 


0.95 


0.7975 


0.95 


0.80 


0.95 


0.80 


0.95 


0.80 


0.95 


0.80 


0.955 


0.80 


0.96 


0.80 


0.96 


0.80 


0.96 


0.80 


0.960 


0.80 


0.97 


0.80 


0.97 


0.80 


0.967 


0.80 


0.966 


0.80 


0.98 


0.80 


0.98 


0.80 


0.974 


0.80 


0.972 


0.80 


0.99 


0.80 


0.99 


0.80 


0.981 


0.80 


0.978 


0.80 


1.00 


0.80 


1.00 


0.80 


0.988 


0.80 


0.984 


0.80 


1.01 


0.80 


1.01 


0.80 


0.995 


0.80 


0.991 


0.80 


1.02 


0.80 


1.02 


0.80 


1.002 






0.80 


1.03 


0.80 


1.03 


0.80 


1.009 






0.80 


1.04 


0.80 


1.04 


0.80 


1.016 






0.80 


1.05 


0.80 


1.05 


0.80 


1.023 







N a 


= 12 


Ns 


= 16 


N s 


= 20 


30000 


30000 


30000 


A 


0A 


A 


0A 


A 


0A 


0.80 


0.865 


0.80 


0.865 


0.80 


0.865 


0.80 


0.870 


0.80 


0.870 


0.80 


0.870 


0.80 


0.875 


0.80 


0.875 


0.80 


0.875 


0.80 


0.880 


0.80 


0.880 


0.80 


0.880 


0.80 


0.890 


0.80 


0.890 


0.80 


0.890 


0.80 


0.900 


0.80 


0.900 


0.80 


0.900 


0.80 


0.910 


0.80 


0.910 


0.80 


0.910 


0.80 


0.920 


0.80 


0.920 


0.80 


0.920 


0.80 


0.930 


0.80 


0.930 


0.80 


0.930 


0.80 


0.935 


0.80 


0.935 


0.80 


0.935 



TABLE III: Parameters sets for PT runs in Fig. [9] left branch. 



for each ensemble for the paths drawn in Fig. [T] are listed 
in the columns of TableslIlto lVIl Remember that the path 
along which we are passing through the finite tempera- 
ture transition at fixed A = 0.8 starts with a horizontal 
piece at fixed (3 A = 0.95. The factor 2 for N in Tables U 
and [IT] refers to the runs with and without modified ac- 
tion S , respectively. In order to remain on the safe 
side in some results we have omitted the first and last 
elements of the ensembles, since the latter have no fur- 
ther configuration to swap with. The respective errors of 
these points might not be of a comparable quality. In the 
literature one can find that by adjusting the parameter 
spacing such that the endpoints get visited with the same 
probability as the neighboring points the errors tend to 
become comparable. 

For the fixed A path of Fig. [fl along which the main 
simulations have been performed, we can fit very well the 



TABLE IV: Parameters sets for PT runs in Fig. |U right 
branch. 
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A 
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0.92 
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0.932 
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0.96 


0.80 
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0.80 
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0.80 
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0.80 
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0.80 
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0.80 
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0.80 
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0.80 


1.04 


0.80 


1.016 


0.80 


1.05 


0.80 


1.05 


0.80 


1.023 



TABLE V: Parameter sets for PT runs in Figs. IFH and [TTI 

step 6(3 a needed to keep u> fixed with a scaling law of 
the form 

6(3 A (cu,N s )^^, (9) 

where we find a(12%) = 2.15(3) in the (3 A = 0.95 - 
1.09 range considered, although we expect it to change 
with the width and location of the the (X, (3a) window. 
Being u> nothing but the tunneling probability among 
twist sectors, it should be proportional to the probability 
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30000 


30000 
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0A 


A 


8a 


0.870 


0.40 


0.830 


0.65 


0.883 


0.40 


0.850 


0.65 


0.888 


0.40 


0.865 


0.65 


0.895 


0.40 


0.888 


0.65 


0.905 


0.40 


0.895 


0.65 


0.915 


0.40 


0.910 


0.65 


0.925 


0.40 


0.925 


0.65 



P A =1 ,X=0.768 



TABLE VI: Parameter sets for two of the parallel tempering 
runs in Fig. [TUJ 



to create a vortex. Since the cost to generate the latter 
should scale with an area law @, the N% dependance 
of the former is easily understood. Eq. implies that 
to explore a fixed region A[3a of parameter space the 
number of ensembles will scale like 



K 



(10) 



As an illustration of the ergodicity of the algorithm 
we show in Fig. [2] the MC time histories of the twist 
observables Zi t i and of the adjoint Polyakov loop for two 
ensembles belonging to the same test PT simulation with 
K = 10 and N s = 16, one below (phase I) and one above 
the bulk transition (phase II). We have first let the sin- 
gle ensembles evolve separately with standard Metropo- 
lis updates, i.e. not updating the permutation table B. 
The barriers among sectors are huge @ and practically 
impossible to overcome without PT within phase II: the 
system is simply stuck in a fixed topological sector. In- 
deed, as Fig. [5] shows the twist variables remained stable 
over 6000 sweeps until we have turned on the full PT 
updates. The system started then frequently to tunnel 
among all sectors. Actually below the bulk transition 
there is no substantial difference between the two algo- 
rithms, since the disorder induced by the Z2 monopolcs 
lets any algorithm be ergodic for the simple reason that 
topological sectors are ill-defined, all twist values fluctu- 
ating around zero. The difference is however dramatic 
above the bulk transition in phase II, where tunnelling 
among well-defined topological sectors is enabled by the 
PT algorithm. 



C. Cross and autocorrelations 

In PT the K ensembles are generated in a correlated 
way. Therefore, the full non-diagonal covariance matrix 
for the observables has to be taken into account. The 
latter is obtained from the general correlation functions 
which, for an observable O and a number N of updates, 
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FIG. 2: MC histories of twist variables z;,4, i = 1,2,3 and 
of the adjoint Polyakov loop La for two PT ensembles. The 
lattice size is 16 s x 4. 



are defined as 



s=l 

N N 

i5>M)(ffS><' 



(ii) 



s" = l 



For j — k they are the usual autocorrelation functions, 
while for j k they describe cross correlations between 
different ensembles. 

The covariance matrix is obtained [52] by using the 
general correlation function of Eq. (jllj) and generalizing 
the derivation in Ref. [H| for the case j = k, which gives 



N-l 



^T,{ 1 -^){ R Mt) + R kJ (t))- (12) 
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The diagonal elements of Eq. (fP2|) are the variances of derivations given in Ref. [53[ for the diagonal case to cal- 



Oj usually written as 



var(0,) = 2r . 



(13) 



culate covariances of covariances from the Rjk(t) data 
only. However, in practice one can hardly get enough 
statistics for this. 



with the integrated autocorrelation times Tj 



D. Correlation results 



When evaluating tj according to Eq. (TTJ} in practical 
simulations the summation up to N — 1 makes no sense 
since pj (t) is buried in the Monte Carlo noise already for 
relatively small t. Therefore, it has been proposed (53l.l54j| 
to sum up only to some smaller value M of t. However, in 
practice such procedure is not stable against the choice of 
M and neglecting the rest is a bad approximation. The 
proposal to estimate the remainder by an extrapolation 
based on the t- values M and M—l [55| is still inaccurate 
in general. A more satisfying procedure is to describe 
the rest by a fit function based on the (reliable) terms of 
Eq. (TT2")) for t < M and on the general knowledge about 
the Markov spectrum. This procedure has led to very 
good results results in other applications [561 ]. 

In order to apply the latter strategy to determine the 
off-diagonal entries in Eq. (|T2)) one has to study how spec- 
tral properties enter the problem. This is possible intro- 
ducing an appropriate Hilbert space [HI, IE?} • Working 
this out, in Ref. [53 ] the general representation 

Rjk(t) = ^ajkrll with |7r| < 1 (15) 

r>l 

has been obtained, where only the coefficients djkr de- 
pend on the particular pair of observables while the eigen- 
values 7 r are universal and characteristic for the simu- 
lation algorithm. To explain the behaviour of the off- 
diagonal elements the approximate functional form 

ii. feW « (£r>lW for |j-fc|<* 1 (16) 
3 W \ for 0<t< \j-k\ J V ' 

has been derived [52| for j ^ k, which indicates a maxi- 
mum at t — \j — fc|. 

For the numerical evaluation of Eq. p2[) the method 
mentioned above is to be used, generalizing it to the off- 
diagonal elements. The fits in the noisy region can exploit 
the universality of the Markov spectrum and the fact 
that after some time only the slowest mode survives. Of 
course, such evaluation is limited by the available statis- 
tics. To calculate errors one has to account for the cross 
correlations between the ensembles. To be able to do 
this one has to rely on fits to the data. The respective 
fit method is well know from the treatment of indirect 
measurements (see e.g. Ref. [58]). For the application to 
PT the details have been worked out in Ref. [52|. To 
obtain errors for the covariances one can generalize the 



Typical examples of correlation functions Rjk{t) (nor- 
malized to Rjj(0)) are shown in Fig. [31 0] and [5] for the 
twist variable zx,4 as introduced in Section 2. Along 
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FIG. 3: Example of a normalized diagonal correlation func- 
tion for 21,4 at /3a = 0.98 and A = 0.8 on a 16 3 x 4 lattice. 
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FIG. 4: Example of a first off-diagonal normalized correlation 
function for 21,4 between ensembles taken at [5a = 0.98 and 
0.99 for A = 0.8 and lattice size 16 3 x 4. 

with all 2^4 we have also used the autocorrelations for 
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FIG. 5: Example of a second off-diagonal normalized correla- 
tion function for 2:1,4 between ensembles taken at (3a = 0.98 
and 1.00 for A = 0.8 and lattice size on 16 3 x 4. 



the Polyakov loop as an additional source to determine 
the eigenvalues "f r in our analysis. For the off-diagonal 
elements Rjk (t) which show a clear signal above the noise 
we generally observe a maximum at roughly t — \j — k\, 
thus seeing indeed the behaviour predicted by Eq. (fl"6|) 
(within errors) in our data. It is usually difficult to iden- 
tify more than two or three off-diagonals above the noise. 
The correlations tend moreover to decrease with increas- 
ing volume. From our data we can clearly conclude that 
the off-diagonal elements of the general correlation func- 
tions are decreasing with the distance from the diagonal, 
their contributions being reasonably smaller than the di- 
agonal one, indicating that cross correlations do not play 
an essential role. 

Fig. [5] shows the integrated autocorrelation times for 
all twist variables 2^4 obtained at different volumes for 
each parallel configuration along the paths at A = 0.8 in 
Fig. [TJ Autocorrelations clearly decrease with the vol- 
ume, as expected. As we shall see, the two peaks corre- 
spond to the bulk transition and the finite-temperature 
transition. 

The obscrvables that will suffer most from correlations 
in PT are obviously those whose expectation value de- 
pends significantly on twist sectors. In our practical case 
only the twists themselves and the vortex free energy are 
sensitive to the ergodicity properties of the algorithm. 
For such observables errors will be given by a combi- 
nation of statistical errors, estimated by bootstrapped 
sampling, and auto/cross-correlation errors given by er- 
ror propagation of the errors on z/j. Other observables, 
like the Polyakov loop and the Pisa disorder parameter, 
are roughly twist independent away from the deep de- 
confined phase, which we anyway do not reach in our 
simulations [III [3l[ . For them only the statistical errors 
will be relevant. 
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FIG. 6: Integrated autocorrelation times for the twist vari- 
able Zi4 at different volumes. The numbers at the horizontal 
axis enumerate the parameter points (A, /3a) for the PT en- 
sembles corresponding to Tab. IIVI 



IV. RESULTS 
A. Monopole condensation 

The computation of the Pisa order parameter p can 
be extended to the parallel tempering approach in a 
straightforward way. Fig. [Jj shows p — log(/i) for 
fixed A = 0.8. As discussed in Ref. [3l|, to prove confine- 
ment through monopole condensation p should be small 
and bounded from below in the confined phase, display a 
dip at the deconfinement transition and reach a negative 
plateau whose value should scale like — 0(N S log N s ) in 
the deep deconfined phase. 

The dip in Fig. [JJ shows the position of the finite tem- 
perature transition. The region left to it, where p should 
roughly vanish, is too close to the bulk transition to ap- 
proach its Pa ~ * value, p indeed has a dip at the 
bulk transition, as shown in Fig. [8] [26|, [3l|, where Z2 
monopoles disappear. Both phases left and right of the 
bulk transition are however still confining as long as p re- 
mains bounded from below on both sides [26|. Fig.[5]com- 
pares the occurrence of the second p-dip around (3a — 0.9 
with the existence of a corresponding peak in the suscep- 
tibility of the average twist as defined in [28| . While the 
location of the latter peak is temperature- independent, 
its height cannot be used to get the scaling with the 4-d 
volume since this should be calculated at T = 0. For the 
susceptibility the latter can be done, obtaining critical ex- 
ponents in accordance with Ising 4-d, as in [281 ]. On the 
other hand, the Pisa disorder operator definition we use 
makes only sense at T ^ (for a definition at T = see 
[HI). Some caution is therefore necessary in interpreting 
the results of Fig. [JJ The region we investigate is, by the 
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FIG. 7: p vs. (5a for A = 0.8 and various lattice sizes. 



(p)erg — 



J SO(3)l 



SO(3)k= 



(17) 



at large /3 a, taking into account the observed 
0(—N s logJVg) plateaus of p for trivial twist and its van- 
ishing at non-trivial twist [H, [3l[ we have (p)erg — 
(/i)| z =o(l — e~^). The latter equation clearly implies 
an exponential vanishing of (p) e rg m the thermodynamic 
limit at high (3 a- 

At low Pa one actually needs a bit more care. In all 
twist sectors p assumes a small constant, bounded from 
below negative value p — > — k ~ —10 [3, HH, therefore 
indicating (p) ^ also for the full ergodic theory, i.e. 
condensation of monopolcs and confinement below (3 A 
[3lj . For every fixed N T (p) e rg can be rescaled post- hoc 
to one through exp(K(3 A (N T )). This rescaling factor will 
necessarily diverge for SU(N), up to logarithmic correc- 
tions, likeiV*" with e = 2(3 (N 2 ~1)/N, [3 = lliV/(487r 2 ) 
being the first coefficient of the /3-function. For SU(2) 
en ~ 1.4. This is of course not an obstacle in normalizing 
(A*) erg = 1 m the limit iV r — *■ oo, although it remains a 
somewhat inelegant feature of the Pisa disorder operator 
in the adjoint formulation. There is however a physical 
motivation for the non vanishing of p, as we will see in 
the following. 



B. Vortex free energy 
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FIG. 8: Comparison between the Pisa order parameter p (up) 
and the susceptibility of the average twist \ z (below) as a 
function of Pa for A = 0.8. 



above exposed limitations of the algorithm, very close to 
both the bulk transition and the physical phase transi- 
tion, so that two competing effects are superimposing. A 
thorough analysis of the whole phase space would be very 
expensive in terms of computer time and could anyway 
hardly be extended to very high (3a- Nevertheless, as ar- 
gued in Ref. [26j], from the fixed twist dynamics of our 
model [H, [3l| one can conclude that for the ergodic the- 
ory the Pisa disorder parameter indicates condensation 
of monopoles in the low (3a region and deconfinement at 
high (3a, provided that a diverging dip at some (3 A exists, 
as Fig. [7] clearly shows. 

Indeed, given that the ergodic expectation value of p 



Having established the physical properties of phase II 
at finite temperature, we will now turn to the 't Hooft 
vortex free energy. As stated above, this observable can 
only be calculated through a fully ergodic simulation. 

In Fig. [5] the free energy of a vortex in lattice units is 
shown as a function of (3a along the A = 0.8 paths of 
Fig. [U The data points start right on top of the bulk 
transition and go on up to slightly above the finite tem- 
perature deconfinement transition. The behaviour up to 
the bulk transition is in agreement with the 't Hooft vor- 
tex argument for confinement: if vortices behave "chaoti- 
cally" (F — 0) then the theory confines (phase I), while as 
deconfinement occurs F ~ oN 2 . As explained above, we 
cannot actually go too deeply into the deconfined phase, 
so we cannot check if for (3a 3> (3 A the data are consistent 
with 0(Ng) plateaus or if they saturate at some value in 
the thermodynamic limit, i.e. if a dual string tension a 
can indeed be measured. To this purpose, assuming that 
the estimate in Eq. (|10p still works at higher (3a, even 
taking into account that for higher volumes the asymp- 
totic behaviour should kick in earlier, we would need to 
simulate around 50 parallel ensembles for each volume, 
again for a statistics of at least O(10 5 ) per configura- 
tion in each ensemble. For volumes with N s > 20, for 
which finite size effects start to be reasonably small, this 
goes beyond the computational power at our disposal, 
although it should be manageable with a medium sized 
PC cluster. A reliable estimate of a would be of extreme 
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FIG. 9: Free energy aN T F v along the A = 0.8 paths of Fig.[T] 

interest in light of the behaviour we find for F in the 
confined phase of phase II, already reported in Ref. [33j . 
Vortex production is there clearly enhanced compared to 
phase I and the free energy stays negative up to the de- 
confinement transition, where it rises to positive values. 

Fig. [TU] shows the free energy in a low (3 a confining 
region well below the finite temperature transition. The 
negative plateau values away from the bulk transition 
at A c ((3a) are consistent with what is observed in Fig. [9] 
and with a vanishing free energy in the limit T — > 
[59l ] . since its value rises again for decreasing (3a after 
reaching a minimum around (3a — 0.65. Larger volumes 
and a better extrapolation would be of course needed to 
confirm this result. 
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FIG. 10: Free energy aN T Fy for varying A — A c at fixed (3a 
and lattice size 16 3 x 4 (see Table Nfy . 

As for the higher twists, although being proper only 



to the toroidal topology they also show a surprising and 
interesting behaviour. Namely, we do not observe a free 
energy proportional to the difference in topological in- 
dex as one would have expected if the twist observables 
Zi t 4 were independent. We indeed observe a strong cor- 
relation among the twist in the different planes, indi- 
cating a non trivial interaction among vortices; as a re- 
sult in the confined phase, although the population of 
the ±1 sectors for the single 2^4 are comparable, the 
distribution of the Vk is such to follow the hierarchy 
^2/3 > u 3 > Vi/3 > u . E.g. for N s = 24, A = 0.8 
and (3 A = 0.96 we find u = 0.020(1), v x = 0.35(5), 
v-2 = 0.48(4) and U3 = 0.15(2). Such hierarchy is quite 
stable with the volume for N s > 16. Taking into account 
that when calculating F v± and i>2 need to be rescaled 
by a factor three to be compared with sectors vq and 
1/3, errors are still to high within the statistics at our dis- 
posal to allow a reliable measure of the free energy for the 
tunnelling other than from/to the sector. Approaching 
and crossing the deconfincment transition the situation 
changes. The trivial sector starts to dominate the par- 
tition function and the free energy to tunnel from one 
sector to another becomes indeed proportional to the dif- 
ference in their topological index. Being the higher sec- 
tors however exponentially suppressed in the deconfincd 
phase their sampling requires longer and longer runs as 
(3a increases; the sampling of sectors with topological in- 
dex higher than one will become in practice eventually 
unfeasible. 



C. C* boundary conditions, monopoles and F 

An interesting alternative check of our surprising nega- 
tive value for F is the evaluation of the vortex free energy 
for the modified action S M needed to define the Pisa dis- 
order parameter in Sec. IIV Al As already discussed in 
Sec. HH C* boundary conditions in the Euclidean time 
direction U(x + aN T i) — U*(x) pose no conceptual prob- 
lem in our adjoint formulation. Any set of twist matrices 
{fi} in the fundamental representation once projected 
onto the adjoint representation become gauge equivalent 
to periodic boundary conditions @. Given that any set 
{£1} is gauge equivalent to the quaternion basis {l2,i&} 
and since C* boundary conditions can be represented 
through the action of = i&2, U*(x) = Q2U (x)Sl\, im- 
posing them makes no difference in the dynamics of twist 
sectors: the configurations that can be assigned to a given 
twist will simply be relabeled with respect to standard 
boundary conditions. In other words the corresponding 
combination of adjoint twist matrices, which would sat- 
isfy a given twist algebra when lifted to SU(2), get reshuf- 
fled by the presence of ia^. A simple listing of combina- 
tions shows however that the number of states leading to 
the assignment of topological sectors z = 0, . . . , 3 remains 
the same. 

C* boundary conditions alone therefore should not 
affect the value of F for S M . There is however the 
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bosonic field giving rise to the abelian monopole through 
its non-trivial transformation property at the boundary 
7r 2 (S77(2)/[/(l)) = tti(£/(1)) = Z (aSHI; its role is far 
deeper and more interesting. For a Lie gauge group G 
with center K and a maximal Cartan subalgebra C it 
is well known that the abelian monopoles classified by 
7r 2 (G/C) = vri(C) will carry center magnetic charges 
classified by ni(G/K) only the ones belonging to the 
kernel of 7Ti(C) — > ni(G/K) will be non singular [60l[6l|. 
For SU(N) this simply means that the N — 1 abelian 
monopoles classified by Z N ~ X will also correspond to 
N — 1 Zjy monopoles; their corresponding Dirac strings 
(sheets in a 3+1 Euclidean formulation) will be in general 
open Zjv vortices. 

In the simple case of SU(2) the assignment is quite 
easy: to an abelian charge n 6 Z will correspond a Z 2 
monopole of charge mod2(w) G Z 2 . Only monopoles of 
charge n = 2k carry no singularity along their world line; 
odd charged monopoles are singular in every gauge. The 
world line of the latter saturates an open Z 2 vortex, while 
only the former are compatible with closed (or no) vor- 
tex sheets [39l |. Our situation is however slightly dif- 
ferent: C* boundary conditions take care of magnetic 
charge conservation by making the unit charge abelian 
monopole being put by hand into the system his own an- 
timonopolc; the Dirac sheet of the Z 2 vortex starting on 
the monopole ends on itself. Such vortex is therefore spa- 
tially closed but cut through by a singularity along the 
world line of the abelian monopole. Moreover such "cut" 
cannot be gauged around; the vortex is "pinned" down 
along one line by the singularity. This should however 
not be detectable by F: the creation of a further vortex 
will bring the system in a state which for our purpose can 
aways be assimilated to either a 2- or 0-twist sector. This 
will make no difference below the bulk transition, where 
the partition function is dominated by open vortices any- 
way and all twist sectors are equivalent: F should still 
vanish there. In the deconfined phase T > T c the cost 
to create a further vortex should still scale like N%, no 
matter what the background is, so that F should become 
again large and positive. Only in the confined phase of 
phase II, if the system prefers the vortex background, 
there should be a difference between S and S : if in the 
former F is negative it should become positive in the lat- 
ter. Fig. QT] shows this to be the case: in the (singular) 
vortex background of S the system prefers to make the 
creation of one further vortex more expensive. Notice 
that the position of the bulk is slightly different for S 
and S* M , e.g. at A = 0.8 it shifts from [3 A ~ 0.89 - 0.90 
to (3a — 0.92 — 0.93; starting our PT ensemble right on 
top of it we are able to span the range of couplings in one 
single run, reaching however the limits of our acceptance 
rate criterion for N s = 20. Why F turns out roughly to 
be minus the half of what it was for S is not obvious. 
It might be related to the non-standard character of the 
S M vortex background or to the effective superposition 
of 0- and 2-twist sectors. In light of a recent proposal 
16211 such observable could also be linked to the direct 
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FIG. 11: aN T F v for S M along the A = 0.8 paths of Fig. Q] 

evaluation of the monopole mass. 

The previous discussion offers also an elegant explana- 
tion both to the sensitivity of p to the bulk transition 
and to its non-vanishing in the low (3a limit of phase II 
observed in [26|, [3l[ and already discussed in Sec. IIV Al 
Action S, consistently with what is expected in the con- 
tinuum limit, completely suppresses the presence of Z 2 
monopoles, allowing only topological, i.e. closed, vor- 
tices. The construction of p we have adopted, following 
[3lL [35l] . uses a unit charged abelian monopole, which as 
discussed above introduces a singularity exactly through 
the presence of a Z 2 monopole. The physical states de- 
scribed by S can therefore never be equivalent to those of 
S M and hence the discrepancy. This was already "pre- 
dicted" in Ref. [39j , where an alternative construction 
using even charged monopoles was proposed as the only 
one fully consistent with the continuum pure Yang-Mills 
action. Singular gauge configurations will not be allowed 
there and only genuinely closed vortices will exist. Such 
modified construction should therefore be the obvious 
choice for adjoint actions and its feasibility might be 
worth to explore. The discrepancy in p observed here 
cannot be detected through actions transforming with 
the fundamental representation: although they should 
disappear in the continuum limit, Z 2 monopoles are still 
present both with and without monopole background 
in the range of parameters commonly used in simula- 
tions, while the topological structure of vortices is any- 
way blurred by the fixed boundary conditions. 

V. ELECTRIC FLUX AND VANISHING OF F 

A few comments are in order to clarify some properties 
of F discussed in the literature which could seem to be 
in conflict with our result. In Ref. [8j an exact duality 
relation between electric and magnetic fluxes implying 
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the vanishing of F was proved under two main assump- 
tions: first one must be able to define a set of regularized 
operators in the fundamental representation; second the 
limit T — > has to be taken after the finite temperature 
compactification of the time direction. As for the lat- 
ter, we actually agree with F = —TAS — ► as T — » 0, 
which is almost obvious in our formulation. The former 
condition has more far reaching consequences. It was ex- 
tended to finite temperature and is indeed the key point 
of the derivations of many interesting dynamical relations 
in Yang-Mills theories, mostly using reflection positivity 
for the fundamental Wilson action on the lattice [63j |. 
In an adjoint discretization however observables like the 
fundamental Polyakov loop, the fundamental trace of the 
electric flux free energy etc. needed to prove such rela- 
tions are all undefined. Formally their expectation values 
and all their correlators vanish identically. 

This can be better understood within the Hamilto- 
nian formulation of lattice Yang-Mills theories, adapt- 
ing the exact construction of the SU(2) Hilbert space 
given in Ref. [641 ]. Since the irreducible representations 
of SO (3) are simply the integer representations of SU(2), 
the Hilbert space of 50(3)-invariant states in each twist 
sector is given by the subset of the SU (2) states described 
in Ref. [6J] (with corresponding t.b.c.) with all links and 
intertwiners labelled by integer spins only. Since topo- 
logical sectors are automatically accounted for in 5*0(3) 



with p.b.c. we have 7iso(3) = ©i^so(3)' * being the 
winding number corresponding to the twist sectors, with 
^SO(3) ^S(7(2) ^ ne Axed twist spaces mentioned above. 
It is now straightforward to see how fundamental op- 
erators annihilate the Hilbert space of SO(3)-invariant 
states. Let us illustrate this by an example. Take a fun- 
damental Wilson loop, winding or not around the bound- 
aries: on the SU(2) Hilbert space such operator can be 
represented as a closed string of spin 1/2 located on the 
corresponding links; as proven in Ref. [64] its action on a 
generic SU(2) state will generate states where the repre- 
sentations labelling the links and intertwiners of interest 
are simply composed with the 1/2 representation via the 
common spin composition rules. Since a state in the 
50(3)-invariant subspace defined above will only carry 
integer spins, its composition with our operator will gen- 
erate states which will necessarily carry semi-integer rep- 
resentations on the links of interest, i.e. it will not belong 
to H.so(3) anymore. Generalizing, the action of a funda- 
mental operator Jon a state tp = J2% •^iV'ti S ^so(3) 
will generate states living in the orthogonal complements 



J-ipi s ^5o( 3 ) with 



0, which is by definition the only state in the 
i,± 



n 



SU(2) 



SO(3) 



SO(3) 



(18) 



tation 1 , as it is the case for Eq. (JTJ), it is legitimate to 
restrict the whole treatment to Hso(S)' Then obviously 

intersection H l so{3) f| Kjo(3)- 

A trivial consequence of the above arguments is that 
with an adjoint action reflection positivity constraints 
can only be invoked for adjoint observables, ensuring e.g. 
that adjoint Polyakov loop correlators will be positive 
definite; constraints derived from fundamental operators 
will be invalid. This concerns e.g. the reflection positi- 
vity constraints among fundamental Wilson loops, static 
quark potential, electric flux free energy and vortex free 
energy (see e.g. Appendix I in Ref. [63]). In particular 
this latter constraint is interesting since it could seem 
to contradict our result for F. Such relation between 
the Fourier transform of the vortex free energy and the 
electric flux free energy, essential to derive the vanishing 
of F in the confined phase, is ill defined in an adjoint 
theory, since it needs the action of a fundamental maxi- 
mal Wilson loop winding around the space boundary to 
be established, as in Eq. (4.6) of [8]. With an adjoint 
weight in the partition function Eqs. (4.9) and (4.10) of 
the above reference become identically zero. Therefore 
the operator there given in Eq. (5.2), although still well 
defined, cannot be related to the projector onto a state of 
fixed electric flux. Also the alternative definition given 
in Ref. [32J through fundamental Polyakov loops hp 
modified via a twist eater at the time boundary is only 
valid in a fundamental theory (with twisted boundary 
conditions). In an adjoint theory the first line of their 
Eq. (15) cannot be inverted. Alternatively, correlators 
of Lp will in general vanish identically giving no useful 
bound on their sign via reflection positivity for the only 
potentially non- vanishing case of maximal displacement. 
There is therefore no guarantee that the right hand side 
of their Eq. (16) will be positive, i.e. that it can indeed 
be interpreted as the exponential of a free energy. 

An electric flux operator can of course be defined also 
in our adjoint model: it will simply be given by an adjoint 
maximal Wilson loop winding around the space bound- 
ary. Such operator cannot however be related to the vor- 
tex free energy as defined in Eq. Q. 



VI. CONCLUSIONS 

In this paper we have studied on the lattice at finite 
temperature a pure SO(3) gauge theory, which trans- 
forms under the actual gauge symmetry group of pure 
SU(2) Yang-Mills in the continuum. Extending the anal- 
ysis in Ref. [33[ we have employed the Pisa disorder op- 
erator for monopolc condensation to establish the prop- 
erties of the theory within the weak coupling phase II, 



When the dynamics can be described by a Hamiltonian 
H = E c + V transforming under the adjoint represen- 



1 Only the transformation properties of V need to be specified; E c 
is always diagonal with our choice of basis |04| . 
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which allows a well defined continuum limit, finding con- 
fined and deconfincd phases separated by a transition of 
presumably second order consistent with the universality 
class of Ising 3-d. 

The vortex free energy F is however found not to van- 
ish at T ^ in the confined phase of SO (3). As discussed 
above, standard arguments for its vanishing up to T c can- 
not be applied when working with an adjoint action, so 
that our result is in itself not contradictory. Moreover, 
the vanishing of F in the confined phase is in general just 
a sufficient condition for confinement. The full adjoint 
theory discussed here possesses neither center symmetry 
nor well defined fundamental observables. Only adjoint 
observables make any sense, a fundamental string tension 
being impossible to calculate. There seems therefore no 
compelling physical reason for the 't Hooft vortex free 
energy to vanish for T ^ 0, since this would be linked to 
an area law behaviour of fundamental Wilson loops, i.e. 
to the existence of a fundamental string tension, which 
however loses its meaning as soon as the specific proper- 
ties of semi-integer discretizations are lifted. Indeed, it 
would be surprising to establish the existence of an order 
parameter for the breaking of a symmetry the theory does 
not possess 65[. This however does not in our opinion 



violate universality, since such properties are not essen- 
tial to describe the dynamics of continuum pure Yang- 
Mills theories. Only physical properties like deconfine- 
ment temperature and universality class or the glueball 
spectrum are preserved independently of the discretiza- 
tion used. Only those should therefore be addressed in 
trying to establish a topological mechanism of confine- 
ment valid both in pure Yang-Mills theories and in full 
QCD. If no symmetry breaking and therefore no order 
parameter is available the properties of the phase transi- 
tion can still be established throu gh t he specific heat or 
other thermodynamic observables jl6l . [l7j . 

The non-vanishing of F implies that the dual string 



tension cannot serve as an order parameter for the adjoint 
theory. For this purpose it would need to vanish exponen- 
tially in the thermodynamic limit in the confined phase, 
while our simulations indicate that it will vanish at most 
as 0(N~ 2 ). Our results show therefore that between 
the monopole condensation parameter and the 't Hooft 
vortex free energy (i.e. a dual string tension) only the 
former seems to retain all its order-parameter character 
independently of the discretization used, at least once 
rescaled to assume a constant value in the confined phase. 
This is a simple consequence of its property to vanish ex- 
ponentially above the critical temperature, as e.g. the 
magnetization for a ferromagnet. The rescaling should 
anyway become unnecessary if one adopts the alterna- 
tive prescriptions given in [39(. Monopole condensation 
can therefore play the role of an order parameter both 
for the SU (2) invariant quenched theory, possessing cen- 
ter symmetry but no vortex topological sectors, and for 
the SO (3) invariant pure Yang-Mills theory, where cen- 
ter symmetry is absent. Vortex topology does not provide 
there a suitable order parameter. 

Why the 0- twist sector gets strongly suppressed in the 
confined phase causing vortex free energy to take a neg- 
ative value, i.e. why the Yang-Mills action finds it en- 
ergetically favorable to create at least one vortex, and 
why twist observables in different directions hint at a 
non trivial interaction pattern among vortices are ques- 
tions which need to be answered and that we will try to 
address in the near future. 
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